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We present here the resuhs for Higgs and Supersymmetry prospective 
searches at the Large Hadron Colhder. We show that for one year at 
high luminosity, Standard Model and MSSM Higgs should be observed 
within the theoretically expected mass range. MSUGRA and restricted 
phenomenological MSSM searches lead to discovery of up to 2.5 TeV 
squarks and gluinos. 
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1 Higgs searches 



1.1 Standard Model Higgs 

The search for the Higgs boson is one of the major task of the LHC Several decay 
channels of Higgs have been explored at the LHC among them : ZZ,WW and 77 
provide the best discovery possibilities. We can see in figure 1 the expected branching 
ratio as function of the Higgs mass. We can notice that bb and WW are the dominant 
mode. 
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Fig. 1: Branching ratios for major Higgs decay channels as a function of its mass. 




Fig. 2: i/ ^ 77 with CMS 
One of the most important and challenging mode is if —> 77 in the Higgs mass 
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range 100 GeV to 150 GeV. It requires excellent calorimetric precision both in en- 
ergy and angular measurement. True 77 production as well as reducible 7jet ( with 
7r° faking photon) should constitute the main background. The figure 2 illustrate 
the diphoton invariant mass as expected with CMS detector for one year at high 
luminosity. 

We can see in figure 3 the expected significance {S/ \fB) for one year of high 
luminosity in ATLAS [1] experiment over the full Higgs mass range. Combining all 
possible channels should give more than 5 cr, even in the low luminosity regime. 
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Fig. 3: Higgs boson discovery with ATLAS, high luminosity 



1.2 Supersymmetry Higgses 

In the Minimal Supersymmetric Standard Model, we expect 5 Higgses, two charged 
and 3 neutral Higgses. For the lightest Higgs, h ^ bb remains the most important 
mode while for heavier Higgs : H/A TT,tt,^ii are interesting modes. As can be 
seen from figure 4 in the case of minimal mixing, the full range m^i = 50 — 500 GeV 
,tanf3 = 1 — 50 can be covered. The most difficult region correspond to intermediate 
rriA and tan(3. In the case of light SUSY particles we would also observe decay of 
Higgs into a pair of sparticles. 
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Fig. 4: MSSM Higgs boson discovery with ATLAS, high luminosity 

2 Supersymmetry searches 

2.1 Introduction 

One of the aims of the LHC collider is to search for the physics beyond the 
Standard Model (SM). One of the direction of this search is a possible discovery 
of superpartners of ordinary particles as expected in Supersymmetric extensions of 
SM (SUSY). SUSY, if it exists, is expected to reveal itself at LHC via excess of 
(multilepton +) multijet + E^*^'' final states compared to Standard Model (SM) 
expectations . 

The main goal of the LHC studies is to evaluate the potential of the CMS and 
ATLAS detectors [2,1], to find evidence for SUSY. It deals, first with a semi- inclusive 
search mainly squarks and gluino mass reach, as the production cross section of 
these strongly interacting sparticles (pair production or in association with charginos 
and neutralinos) dominates the total SUSY cross section over a wide region of the 
parameter space. These studies were done first in the framework of MSUGRA model 
and at a later stage in a constrained MSSM. 

Finally we will consider more exotic models like R-parity violated and Gauge 
Mediated Supersymmmetry Breaking models. 
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2.2 MSUGRA Studies 



The large number of SUSY parameters even in the framework of Minimal extension 
of the SM (MSSM) makes it difficult to evaluate the general reach. So, for this study 
we restrict ourselves at present to the mSUGRA-MSSM model. This model evolves 
from MSSM, using Grand Unification Theory (GUT) assumptions (see more details 
in e.g. [3]). In fact, it is a representative model, especially in case of inclusive studies 
and reach limits are expressed in terms of squark and gluino mass values. 

The mSUGRA model contains only five free parameters : 

- a common gaugino mass (mi/2) ; 

- a common scalar mass (mo) ; 

- a common trilinear interaction amongst the scalars (Aq) ; 

- the ratio of the vacuum expectation values of the Higgs fields ; 

- a Higgsino mixing parameter n which enters only through its sign {sign{ii)). 

2.2.1 Simulation procedure 

The PYTHIA 5.7 generator [4] is used to generate all SM background processes, 
whereas ISA JET 7.32 [5] is used for mSUGRA signal simulations. The CMS JET 
(version 4.51) fast MC package [7] is used to model the CMS detector response, since 
it still looks impossible to perform a fuU-GEANT [6] simulation for the present study, 
requiring to process multi-million samples of signal and SM background events. 

The SM background processes considered are : QCD 2^2 (including bb), tt, 
W + jets, Z + jets. The pT range of all the background processes is subdivided into 
several intervals to facihtate accumulation of statistics in the high-^r range : 100-200 
GeV, 200-400 GeV, 400-800 GeV and > 800 GeV (additional interval of 800-1200 GeV 
is reserved for QCD). The accumulated SM background statistics for all background 
channels amounts to about 200 millions events. 

The grid of probed mo, mi/2 mSUGRA points has a cell size of Amo=Ami/2=100 
GeV for mo <1000 GeV and Amo=200 GeV, Ami/2=100 GeV for mo >1000 GeV. 
This was also probed with the appropriate mixture of signal and pile-up events. The 
kinematics of signal events is usually harder than that of SM background for the 
interesting regions of maximal reach of squark-gluino masses . The cross section of 
the background is however higher by orders of magnitude and high-py tails of different 
backgrounds can have a kinematics similar to that of the signal. 

In figure 5 we compare some signal distributions for the point (mo, mi/2) = 
(1000,800), corresponding to m^ m^^ 1900 GeV, m^^o = 351 GeV, m^o = m-± = 

668 GeV, and distributions of the sum of all SM background processes for the E^J*''* 
signature. 
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Both signal and background histograms contain only events satisfying first level 
selection criteria {Nj^ts > 2 and i?™*'''^ > 200GeV). Only the hardest jet and lepton 
in the event are shown in distributions in figure 5. 




Fig. 5: Comparison of the mSUGRA signal and SM background, with 100 fb , in one point 
of mSUGRA parameter space : mo=1000 GeV, mi/2=800 GeV (m^ssm^-^ps 1900 GeV) for 
the E™'^'' signature. Initial cuts are applied. 

Since the topology of signal and background events is rather similar already after 
first level selection cuts, the difference in the angular distributions and circularity is 
not significant either, it is thus not very useful to apply cuts on these variables too. 
The difference in the lepton p^ distributions is also not very pronounced as signal 
leptons are produced in cascade decays, thus loosing "memory" about the hardness of 
the original process. But for extremely high masses of squarks or gluinos (~ 2 TeV), 
there is some difference in the angular and p^^ distributions between signal and the 
total SM background. One can deduce from figures 5 that cuts on the jet multiplicity 

and E^*'*'' are the most profitable ones for background suppression . Of course, 
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there is inevitable correlation between variables both in signal and background, e.g. an 
obvious correlation between E^*** and the hardest jet in QCD events, since there 
-^mtss niainly produced by neutrinos from b-jets and/or high-E:j^^* mis- measurement. 
This can lead to a degradation of the efficiency of some cuts, if fixed cuts are used. It 
is thus more profitable to have adjustable cuts to meet various kinematical conditions 
in various domains of mSUGRA parameter space and take into account difference in 
topology between various signatures. For this reason, we search among 10000 cuts 
combination, the best selection at each generated MSUGRA point. 

2.2.2 Results 

Figures 6 and 7 show the main results of our study for mSUGRA assuming an 
integrated luminosity of 100 fb~^. Figure 6 contains isomass contours for squarks (q), 
gluino (g) and lightest scalar Higgs (h). Numbers in parenthesis denote mass values 
of corresponding isomass contour. The neutralino relic density contours from ref. [8], 
for mSUGRA domain mo<1000 GeV, mi/2<1000 GeV, are shown for ilh^ = 0.15, 
0.4 and 1.0. Value Q/i^>l would lead to a Universe age less than 10 billion years old, 
in contradiction with estimated age of the oldest stars. The region in between 0.15 
and 0.4 is favoured by the Mixed Dark Matter (MDM) cosmological models. It is 
a rather general situation that for all investigated sets of mSUGRA parameters the 
best reach can be obtained with the E^*''* signature. The more leptons required - 
the smaller reach, as can be seen from figures 7. The cosmologically preferred region 
flh^<0.4 seems to be entirely within the reach of CMS. In figure 6 we also show our 
calculations for the E™"*' signature reach for an integrated luminosity of 300 fb~^, 
trying to estimate the ultimate CMS reach. 
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Fig. 6: 5 sigma reach contours for various final states with 100 fb~^ (see also comments in 
text). 
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Fig. 7: Simplified figures gathered together. 
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2.3 MSSM 



2.3.1 A restricted pMSSM 

The model used in the next study is a phenomenological MSSM (19 parameters), 
but with a further reduction of the number of the parameters which allows us to 
perform simpler simulation, while keeping the diversity of the signatures of MSSM 
events . We take into account, respectively, the mass unification of squarks and slep- 
tons (universality of the three particles generations) and the unification of trilinear 
couphng. This leads to 9 free parameters : tan/3. Ma, A*, Mi, M2,M~, M~, Mj, As- 

2.3.2 Signal 

The goal of this study is to evaluate the capacity of CMS detector to observe 
MSSM signals. Two reasons urge us to study this pMSSM after the MSUGRA study. 
First , MSUGRA is a rather constrained model with only five free parameters. Con- 
trary to MSUGRA, the pMSSM has no fixed hierarchy of masses as shown later. 
Moreover, we try to estimate MSSM parameter values using kinematical quantities 
measured by CMS, using a fast simulation program. We use a model with 9 param- 
eters, which constitute an hyperspace with 9 dimensions. In order to facilitate the 
analysis, we use a discretization of the parameters. The choice of the number of value 
of each parameters depends of the parameter sensitivity. We used a grid for squarks 
and gluinos masses with 9 values evenly spaced between 600 and 3000 GeV, because 
the events characteristics at LHC depend primarily on these two masses. On the other 
hand, the observables are not very dependant on the parameter tan j3 and thus we 
only use two values in order to distinguish the larges ones and the small values from 
it parameters. Thus, the values which we selected for each parameter of this analysis 
are : 

- Mj: 200, 1000, 3000 GeV 

- Ml : 100, 500, 1000, 2000Ge\/ 

- M2 : 100, 500, 1000, 2000Ge1/ 

- Ma : 200, 1000, 3000Ge1/ 

- A3 : 0, 2000 GeV 

- fi : 200, 500, 2000 GeV 

- tan/3 : 2, 50 

- M~ : 600, 900, 1200, 1500, 1800, 2100, 2400, 2700, 3000 GeV 

- M~ : 600, 900, 1200, 1500, 1800, 2100, 2400, 2700, 3000 GeV 

We end up with a total of 140000 different combination of parameters for each of which 
we generate 1000 events, a compromise between the limits imposed by handling the 
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Fig. 8: Example 1 : distribution of signal and background for differents observables 

data flow and sufficiently small statistical errors. The theoretical and experimental 
constraints make it possible to reduce the number of combination to a total of 17.10^. 
The background production was estimated using Standard Model events leading to 
similar signatures as MSSM events, in order to study the possibility of extracting 
SUSY signals from the background. We illustrate in the following figure the back- 
ground and signal distribution for various variables. 

The values of MSSM parameter for the figure 8 are : 
Mj = 1000 GeV, Ml = 100 GeV, Ma = 500 GeV, M~ = 1800 GeV, M~ = 1800 GeV, 
Ma = 200 GeV, tan /? = 50, /x = 2000 GeV, A3 = 2000 GeV. 

Despite the corresponding low cross section a = 77/6, the significance of this param- 
eter set is, after application of optimisation cuts, equal to 21. 

2.3.3 case of large and close hierarchy of masses. 

For the first example (figure 9) the masses of neutralinos are much lower than 
the masses of squarks, gluinos and sleptons, the production will be dominated by 
neutralinos and charginos. The parameter values are the following : 
Mj = 2000 GeV, Mi = 500 GeV, M2 = 500 GeV, M~ = 2000 GeV, M~ = 2000 GeV, 
Ma = 1000 GeV, tan p = 50, = 200 GeV, A3 = GeV. 

The cross section of this set of parameters is cr = 1.22pb and despite the abundance 
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Fig. 9: Example 2 : distribution of signal and background 



of neutralinos, the low production rate of gluinos and squarks allows nevertheless to 
obtain a significance equal to 10.2. The cuts for figures 9 and 10, are those giving the 
maximum significance (among the list of the optimization cuts). 

For the second example (figure 10) the parameter values are the following : Mj = 
1500 GeV, Ml = 940 GeV, Ma = 2000 GeV, M~ = 1000 GeV, M~ = 1020 GeV, Ma = 
1000 GeV, tan /3 = 50, /i = 1050 GeV, A3 = GeV. The masses of neutralinos, gluinos, 
squarks and sleptons are comparable. The main production proceeds via gluinos and 
squarks with a cross section a = 2.0 pb and a significance equal to 36.3. 

For the third example (figure 11) the parameter values are the following : 
Mj= 1520 GeV, Mi = UbOGeV, M2 = 2000 GeV, M~= 1500 GeV, 
M~ = 1520 GeV, Ma = 1000 GeV, tan /5 = 50, /i = 1500 GeV, A3 = GeV 
The masses of neutralinos, gluinos, squarks and sleptons are comparable but with 
a higher value. The main production proceeds always via gluinos and squarks with, 
in this case, a cross section a = 0.126 pb and a significance equal to 3.2. 

Even for the sets of parameters which would seem difficult (hierarchy of very close 
mass or on the contrary very separeted), this method make it possible to obtain good 
results, but we observe a limitation in the case of close hierarchy of masses with a 
dicovery limit of about 1.5 TeV instead of 2.5 TeV in the other case. 
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Fig. 10: Example 3 : distribution of signal and background 




Fig. 11: Distribution of signal and background 
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Fig. 12: Graph of exclusion in the plane m~vsm~ 
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Fig. 13: limit of discovery of MSSM in plane m~vsm~ 
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2.3.4 Limit of discovery 



Wc now will generalize this study in order to define limit of discovery for the 
MSSM. The significance was calculated for high luminosity / Ldt = 100/6~^. The 
isocurves of significance equal to 5 are given in Figure 12 in the plane m~ versus m~, for 
the set of 7 differents parameters : Mj^ 3000 GeF, Mi = 100 GeV, Ms = 2000 GeF, 
Ma = 200 GeV, tan/5 = 50, /i = 2000 GeK, A3 = OGeV. Each of these isocurves 
represent a specific cut on the number of leptons per event. The curves labelled 0/, 
11 and 21 represent significance equal to 5 in case of each event taken into account 
in the calculation of signal and background possesses respectively 0, 1, 2 lepton(s). 
Figure 13 show the zone of discovery in plane m-vs m~ for all the sets of parameter. 
There is a total of 3750 from the 17000 combinations of parameter which wouldn't 
be discovered. This lead to the limit of discovery of about 2 TeV. 

2.4 Conclusion 

We demonstrated the possibility to discover a phenomenological MSSM using an 
inclusive study in the MSSM parameters space. We can note, at the end of our study, 
that we observed httle difference between MSUGRA and the pMSSM. The hmit of 
discovery correspond to the hmit of the cross section (2.7 TeV at CMS). The only 
difference appears for some points having a specific mass hierarchy. For example in 
the case of close masses, the limit we expect is about 1.5 TeV. 



3 R parity violation 

The conservation of this quantum quantity {R = (— 1)3(-B-l)+25^ allows to conserve 
baryonic and Icptonic number. Nevertheless there is no theoretical reason to impose 
this conservation. We then expect several new terms to appear in the supersymmetric 
lagrangian depending on a A coupling. Experimental constrains lead to very small 
value for the factor (usually less than A < 10~^) multiplying the lagrangian terms 
violating the R parity. Phcnomenologically for (A < 10^^) we expect the LSP to decay 
outside our detectors thus behaving like R-parity conserved signature. We focused here 
on the intermediate situation (10^^ < A < 10^^). The phenomenological behaviour 
correspond to R-parity conserved situation with the decay of the LSP in 3 jets or 
3 leptons or 2 jets and one lepton. Figure 14 illustrates the 3 leptons case where 
the invariant dilepton mass is plotted. As is expected the edge of this distribution 
correspond to the LSP mass. This correspond to MSUGRA model with low mo mass. 

The expected reach for the 3 leptons case in the MSUGRA plane is illustrated in 
the following figure 15. We can note that the mass reach is roughly the same as in the 
R-parity conserved one. In the following figure 16, is plotted the invariant dijet and 
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Fig. 14: Invariant dilepton mass for A121 > 
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Fig. 15: MSUGRA mass reach for one year of low luminosity for A121 or A233 non zero 
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Fig. 16: invariant dilepton and jet mass for A' non zero. 



lepton invariant mass corresponding to A' non zero. We can observe a reconstructed 
peak at 130 GeV matching the LSP mass. 



4 Gauge Mediated Supersymmetry Breaking 

In these model, SUSY breaking has its origin in a separate sector at a relatively 
low scale Fm = (lO^^Ge^)^ (compare to MSUGRA model). An important caracteris- 
tic of this model is to lead to a very light gravitino mass much lower than the other 
spartners. In the minimal GMSB model 6 parameter are considered : Mm messenger 
scale, A = Fm/Mm, A^5 number of messenger family, tan/3, sgrifi, Cgrav ratio of the 
gravitino mass to its value if the only source of SUSY breaking is Fm- For LHC ex- 
periment it is extremely interesting to evaluate the response of the detector to the 
signature provided by these models as they present different phenomenology. We ex- 
pect 2 possible scenarios : 

- X?^G + 7 (X?NLSP) 

- fR^G + l {In NLSP) 
with the following signatures : 
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Fig. 17: Dilepton invariant mass reconstruction for slow ionizing slepton 

In figure 17, we show the channel Xi ^bJ- where we reconstruct a long living 
right slepton using TOF measurement as measured in the muon chamber. We observe 

3 peaks corresponding to Xi)X2 xl- Figure 18 show the channel II — > X?^ ~^ ^R^^ '■ 
peak and — > z// — » Ijiull. Combine previous Xi invariant mass with any of 

4 hardest jets provide the reconstruction of left squark (generated rriq = 648 GeV, 
reconstructed 632 GeV) as can be seen in figure 19. From this study we can deduce 
that the full reconstruction chain can be performed. 

4.1 Conclusions 

The main conclusions of our studies are the following : Standard Model Higgs as 
well MSSM Higgs should be visible over the complete theoretical expected range using 
various decay channels. The SUGRA model, investigated for SUSY search, would be 
detectable through an excess of events over SM expectations up to masses m^ ~ m^ 
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Fig. 18: Trilepton invariant mass 
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Fig. 19: Trilepton and jet invariant mass 
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~ 2.5 TeV with 100 fb^^ This means that the entire plausible domain of EW-SUSY 
parameter space for most probable values of tan/3 can be probed. Furthermore, the 
S/B ratios are > 1 everywhere in the reachable domain of parameter space (with 
the appropriate cuts) thus allowing a study of kinematics of g, g production and 
obtaining information on their masses. The cosmologically interesting region Qh'^ < 
1, and even more the preferred region 0.15 < flh"^ < 0.4, can be entirely probed. 
The mass reach is up to q,g — 2.5TeV and the cosmologically interesting region is 
covered by CMS/ATLAS. Sleptons can be observed up to 400 GeV, LSP up to 400 
GeV. LHC can perform exclusive studies : gluino, neutralino reconstruction and de- 
termine MSUGRA parameters. MSSM investigation show similarities with MSUGRA 
but some specific mass hierarchy might be difficult to explore. LHC can also strongly 
contribute in GMSB model and provide some insights for parameters determination. 
R parity violated ( A < 10~^) studied ( Xi mass determination). These models lead 
to various experimental signature with which LHC detectors can cope. 
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